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SUMMARY 

The axial and. radial fields at any point inside or outside a finite 
solenoid with infinitely thin walls are derived. Solution of the equa- 
tions has been obtained in terms of tabulated complete elliptic integrals. 
* °r the axial field an accurate approximation is given in terms of 
elementary functions. Fields internal and external to the solenoid are 
presented in graphical form for a wide variety of solenoid lengths. 


INTRODUCTION 

The recent great Interest in plasmas either as a source of energy 
or as a propulsion device has resulted in a greatly renewed interest in 
the magnetic fields produced by various configurations of electromagnets 
Of the possible methods of plasma confinement by far the most promising 
appears to be the use of magnetic fields (ref. l). 

The calculation of the fields generated by various electromagnetic 
configurations such as loops, finite helical solenoids, and infinite 
solenoids has been treated by the early classical physicists, but only 
the simplest cases such as the single loop have been calculated for the 
entire field both inside and outside the loop (e.g. , ref. 2). In other 
cases such as the helical solenoid or the finite solenoid the calcula- 
tions have been limited to the axis (ref. 3). Derivations of the off- 
axis positions have been done by Foelsch (ref. 4), and the solutions 
are obtainable by means of a large number of approximate expressions 
which are valid over restricted ranges of size or position. The princi- 
pal difficulty in the calculation of the fields of nearly all config- 
urations has resulted from the fact that the integral solution cannot 
be achieved without the use of various elliptic integrals. Even though 
many of these are tabulated, the calculations involved are laborious. 

Such calculations can, however, be made using modern high-speed com- 
puters since machine programs have or can be written for many of the 
elliptic functions. 
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The purpose of this report is twofold: (l) to derive the equations 

of the axial and radial field at any point within or outside a finite 
solenoid in terms of standard tabulated functions and (2) to plot these 
fields for a number of solenoids. 

SYMBOLS 

A q magnetic vector potential component in ©-direction 

a coil radius 

radial and axial magnetic induction component 
E complete elliptic integral, second kind 

i current in each filament 

K complete elliptic integral, first kind 

k ^/4ar/[| 2 + (a + r) 2 ] 

L coil length 

n number of turns per unit coil length 

r, 0, z cylindrical coordinates 
A Q Heuman lambda function 

p permeability 

z ± 

cp tan 



DERIVATION OF EQUATIONS 


Consider a solenoid as shown in the following sketch: 
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The magnetic field due to this coil is given in terms of the vector 
potential A by 

B = VXA (l) 


where, for the geometry assumed, only the A q component can be nonzero. 
Then equation (l) yields simply 


SA/ 

B r= “ST 


Br. 


x S(rA 0 ) 


( 2 ) 


r Sr 

For a single circular filament, one has 


A 0 


4 it 


_ ui a cos 0 d0 


R 


where R is the distance from the local point on the filament to the 
field point. For a solenoid made up of a series of n filaments per 
unit length, we have then 


An - ^ 

9 2 it 


^*L/2 /• 

I n 

-L/ 2 * / 0 


cos 0 d0 


V(z - Vr^ + a^ - 2ar cos 0 


or 


, _ ajin 

- 2? 


cos 0 d0 


- f + ^ f —= 

Jt Vs + r^ + a^ - 2ar cos 0 


(3) 


where £ = z -7,, ^ = z ±, L/2, and l is the axial distance from the 
origin to the filament. On integrating with respect to £, this becomes 


a = aoni 

A Q 2n 


- J* cos 0 lnj^ + 
0 


+ r^ + &r - 2ar cos 01 d0 (4) 


>]!: 


A more convenient form can be found by integrating by parts: 

|0=« 




a^ni 
2 it 


jsin 9 lnjl + yft 

f 


2 + + a^ - 2ar cos d 


) ^ + 


e=o 



ar sin^0 d0 

^ + V? 2 + r 2 + a 2 

- 2ar cos ©jVl^+r^ + a^- 2ar cos 0 
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The first term vanishes. On multiplying the integrand hy 


4- r 2 + a 2 - 2ar cos g - | 
^2 + r^f ”a^ - 2ar cos 0 - £ 


rearranging terms , and observing that 


use of the limits eliminates one term, there follows 


k G - 


a 2 unir 


ft 
* / 0 [_( 


£ sin 2 © d© 


(a 2 + r 2 - 2ar cos ©)' y£ 2 + r 2 + a 2 - 2ar cos ©l| 


The two magnetic-field, components can now be easily obtained. The radial 
field is found by differentiating equation (3) and yields 


B r ~ " c)z 


/[- 

J o L 1 


cos G diO H 

f^2 + r 2 + a 2 _~2ar cos 0 


To get B z , first evaluate SA^/dr from equation (4). The result, pro- 
ceeding as in the case of obtaining equation (5), is 


^ k Q _ _ auni 

hr 2 * 


a 


^ cos ©(r - a ^os ©)d© 


(r 2 + a 2 - 2ar cos ©) Vl 2 i- r 2 + a 2 - 2ar cos © 


If equations (5) and (7) are put into (2), th2 result is 


B„ = 


^ rA ©) _ auni 


/[■ 

J 0 L 


£(a - r cos g)d© 

( r 2 + a 2 - 2ar cos ©) ^ 2 + r 2 +a 2 _ 2ar cos 0 | 


Equations (6) and (8) describe the magnetic field due to a finite sole- 
noid. Numerical results can readily be found by integrating these equa- 
tions on a computer. However , the results cai also be expressed in 
terms of standard elliptic integrals, which aus already tabulated. This 
we proceed to do. 


Consider B r . This can be evaluated by use of formulas 291.03 and 
412.01 (noting the special case of a 2 = k 2 ) xf reference 5. One has 
successively. 
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where 


k 


2 


4ar 

£ 2 + (a + r ) 2 


( 10 ) 


In a similar manner, B z (eq. (a)) can be reduced to standard elliptic 
integrals. First change the variable of integration to t = cos Q and 
then use formulas 233.19 and 413.06 of reference 5 . There follows* 
successively 


B = -J^ 11 

4jt \/2ar 


f) 


5 


(M 


dt 


(4^ - *) - **>(* 


t l l/(i - 1 


B, 


jani 


2it(a + r) 




5k 


/ 

0 


2 r 9 

l - * sn^u* 

a + r 

4ar 9 

1 - — sn^u 


du 




(a + r ) 2 

^fe KW + Tfe^T A °H s U 


where 


(id 


cp = tan ( 12 ) 

As before, k is given by equation (lO). The Heuman lambda function 
A 0 (cp,kj is tabulated in references 5 and 6 . 


For many purposes, it is convenient simply to know the variation 
of the fields near the axis. As r -» 0, equations ( 9 ) and (ll) reduce 
to the following well-known expressions? 

•n _ [ini | a 2 r 1 

r= 4 L(T 2 + a 2)3/ 2 J| 


(13) 
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B, = ttsi 

z 2 


VP 


: 15 + 

— 

+ a 


( 14 ) 


A convenient approximation for B z , valid whenever r < a and accurate 
to 1 percent in this range, is 


Bz = + 2k' ) ^ + [ p m ' > + 2k’J sin cp| 




(15) 


where m = (l - k 1 )/(l + k 1 )^ k 1 = Vl - Equation (l5) reduces ex- 

actly to equation (14) at the axis. 
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CALCULATIONS 

Equations (9) and (ll) are readily written in dimensionless form 
with the distances given in units of the solenoid radius. Then equa- 
tions (9) to (14) still hold hut with 1, r/a, l^/a replacing a, r, 

^ throughout. 

Plots of the dimensionless axial and radial fields, +4B z /|ani, 
-4B r /(ini, are given in figures 1 and 2, respectively. Calculations were 
made for the ratio of coil length to radius in the range from 1 to 25. 
Note that in figures 1 and 2 the radial distance is given in terms of 
the coil radius (r/a) and, the axial distance is given in terms of the 
coil half-length (2 z/l). 


Discussion 

The figures clearly show that increasing soLenoid length decreases 
the radial variation of the axial field. This result is expected since 
an infinitely long solenoid has a uniform field throughout. For short 
solenoid lengths (fig. 1(a)), the axial field increases rapidly from 
the center to the wall for positions near the center of the solenoid. 

In fact, at the center the curve approaches very closely that for a 
simple loop. 

It should he noted that the radial field is always infinite at 
2z/L =1.0 and r/a = 1. This point corresponds to the edge of the 
current sheet and would he expected to produce such a result. 

In general, uniform fields with total variations of 1 percent can 
he achieved over as much as 60 percent of the internal volume of the 
solenoid if the length is 25 radii or greater. 


-Qnn 
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The calculations presented herein are limited to solenoids with in- 
finitely thin walls but the results can readily be used to find approxi- 
mate solutions for various shaped solenoids of finite thickness with al- 
most any current distribution. Since superposition principles apply, it 
is only necessary to approximate any odd-shaped solenoid by a number of 
thin- walled solenoids and add the fields resulting from each. The accu- 
racy of the answer is, of course, dependent on the number of separate 
solenoidal rings used to approximate the actual shape. 

It is interesting to note that the results obtained herein for mag- 
netic fields are closely related to the velocity fields produced by a 
lifting helicopter rotor (e.g., ref. 7). The physical model is the same 
but the detailed methods of solution are widely different. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 23, 1960 
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Dimensionless axial field, 4B /|ini 
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(b) L/a = 2. 

Figure 1. - Continued. Dimensionless axial field of a finite solenoid. 
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(d) L/a = 4. 
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Dimensionless axial field, 4B Z /Vni 

(e) L/a = 5. 

Figure 1. - Continued. Dimensionless axial field of a finite solenoid. 


006-3 







Dimensionless radius 


E-900 



3.2 


£ .8 1.2 1.6 2.0 2.4 2.8 
Dimensionless axial field, 4B 2 /jj,ni 

(f) L/a = 10. 

Figure 1. - Continued. Dimensionless axial field of a finite solenoid. 
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Figure 2, 


Continued. Dimensionless radial field of a finite 
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(c) L/a = 3. 


Figure 2. 


- Continued. Dimensionless radial fielc of a finite solenoid. 
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Figure 2. - Concluded. Dimensionless radial field of a finite solenoid. 
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